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Mitochondrial DNA Sequences Imply Anhui Musk Deer 
a Valid Species in Genus Moschus 
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Abstract: Anhu: musk deer | Moschus ( moschiferus ‹ berezovkit ) anfuuensis | bas been a taxonomic mystery since its 
discovery in early 80's Ја this paper. with museum samples. we sequenced the complete cytochrome b gene of five An- 
hui musk deer. When compared with other species in Genus Moschus, Anhui musk deer showed a rather level of se- 


quence divergence from all the other species in this genus. The phylogenetic trees constructed by multiple methods sup- 
ported the same topology. in which the monophyly of Anbui musk deer was clearly demonstrated. Therefore. our molecu- 
lar data suggest a full species status tor Anhui musk deer ( Muschus anhuiensis |. rather than a subspecies of either 


М. moschiferus ог М. berezovskii previously suggested by morphological studies. 


Key words: Musk deer: Genus Moschus; Anhui musk deer; Cytochrome b gene; Molecular phylogeny 


Musk deer (Genus Moschus ) are distributed in 
China and ity adjacent areas ( especially the Qinghai- 
Tibetan Plateau and Himalayan areas) (Groves et al., 
1995). Due to deforestation and hunting, musk deer 
populations are decreasing sharply and some species 
are seriously endangered and difficult to find in the 
wild. Taxonomically, musk deer was classified as a 
family ( Moschidae) in pecora group. However. there 
are «till arguments on the number of species and the 
phylogenetic relationships among them (Ellerman & 
Morrison-Scott, 1950; Gao, 1963, 1985; Li, 
1981: Grubb, 1982; Groves & Feng, 1986; Groves 
et af., 1995; Ohtaishi & Gao, 1990; Sheng. 
1989; Wang et al., 1993). Based on the character- 
istics of external and skull morphology. and used in a 
multivariate analysis for measurements, Groves et al. 


(1995) suggested that there are 5 species in musk 
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M . berezovskü . 


M fuscus, M .chrysogaster and M . leucoguster. This 


deer. which are М. moschiferus . 
opinion was confirmed by recent molecular data (Su et 
al.. 1999). 

The so-called “Anhui musk deer" live in Anhui 
Province of China. lt was firstly reported in early 80's 
by Wang et al. (1982) and described as a subspecies 
of M.moschiferus ( M.m.anhuiensis ), and later 
Groves & Feng (1986). classified it as a subspecies of 
M . berezovskii (М. anhutensis ). However. Groves 
et al. (1995) pointed out that; “a skull of Animi 
musk deer falls outside the range of both M . berezovskzi 
and M . moschiferus samples". implying an uncertainty 
ol taxonomic status for Anhui musk deer, In order to 
clarity the ambiguity, we collected the type specimens 
ol Anhui musk deer, extracted DNA from degraded 


leather samples and sequenced complete cytochrome b 
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(Cytb) genes. The molecular data is not only provid- 
ing new taxonomic evidence, but also would he useful 
in cunservation and genetic management of this endan- 


gered group. 
1 Materials and Methods 


1.1 Sample collection 

All the five individual samples were collected trom 
cullections in Department of Biology, Anhui Universi- 
ty. China. They are leather specimen sampled about 
20 years ago. The specimens were morphologically i- 
dentified by analyses of the whole animals. There are 
опе type and one paralype specimens in the five speci- 
mens studied (laheled No.7903 to Nu. 7907 respec- 
tively}. 

1.2 DNA extraction, PCR and DNA sequencing 

For DNA extraction, we followed the method of 
Walsh et ai. (1991) with a few alterations. The de- 
tailed procedures were described in our previous paper 
i^u et al.. 1999}. 

We used a set of primers to amplify Cytb gene of 
mitochondrial DNA (Su et a?., 1999). Each pair of 
primers was designed to cover 200 — 300 bp fragment» 
in Cytb gene. The sequences ul all the fragments were 
then aligned and overlapped to get complete sequences 
of Cytb. PCR was done in the following conditions 
(Robo-Cycler, Stratagene}: predenature at 94°C for 2 
minutes, then cycling at 94% (30 seconds}, 42 - 
50%. (30 seconds), and 72 (I minute) for 40 cy- 
cles, and final extension at 727C for 5 minutes. The 
PCR products were purified bv LMP argrose gel elec- 
trophoresis and roughly quantified by eves through EB 
staining . 

For sequencing, an automatic DNA sequencer 
(ABI Model 377) was used to sequence the double- 
stranded PCR products directly. Cycling sequencing 
method with FS kits (ABI) was employed following the 
instructions of the producer. All amplified fragments 
were sequenced both with the forward and reverse 
primers . 

1.3 DNA sequence alignment and phylogenetic 
analysis 


DNA sequence alignments were done by eyes. 


For phylogenetic analysis, we used multiple algorithms 
for tree construction, including the parsimony method 
in PAUP 3.0 (Swofford, 1989), the maximum likeli- 
hood method in PHYLIP 3.5 ( Felsenstein, 1993), 
and the neighbor-joining method in MEGA 1.02 ( Ku- 
mar et al., 1993). Confidence values for internal 
lineage in parsimony analysis were assessed bv buot- 
strap (Felsenstein, 1985). The Су sequences of 
other species of musk deer were from our previous 
study (Su et al., [999). The published Cytb se- 
quences of cheyrotain ( Tragulus napu ? and moose 
(Ades alces ) were used as outgroup data ( Irwin et 


al.. 1991}. 
2 Results and Discussions 


2.1 Complete Cytb gene sequences and genetic 
divergence 

The complete Cytb gene sequences of the five An- 
hui musk deer were obtained by overlapping the 8 frag- 
mented sequences of each individual. The whole 1 140 
bp sequences of all samples of Ánhui musk are quite 
similar. Onlv two haplotypes were observed with one 
transitional difference (see Fig.1). However, more 
substitutions were observed than expected between An- 
hui musk deer and the other musk deer species, rang- 
ing from 76 to 88 (6.7% – 7.7%) (Table 1). A- 
mong them, 31 substitutions (29 transitions and 2 
transversions) are unique for Anhui musk deer. As 
compared with the divergence level among species and 
subspecies of musk deer, the sequence divergence 
(7.3% on average) between Anhui musk deer and the 
others is apparently heyond the subspecies level. For 
example. the sequence divergence between two sub- 
species of АЎ. beresouskii ( M.b. bijiangensis and 
M .b.caobangis) is only 1.01% (Su et al.. 1999}. 
2.2 Taxonomy and phylogeny 

Using multiple tree-constructing methods ( parsi- 
mony, distance and maximum likelihood), we updat- 
ed our phylogenetic tree of Genus Moschus (Su et al., 
1999). Asa result, all the methods apphed for tree- 
construction support the same topology shown in Fig.2 
with very high hootstrap values. As we mentioned 


above , the genetic divergence between Anhui musk 
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Table 1 Pairwise comparison of seqnence divergence in musk deer species 





species 1 2 3 4 3 б 7 B g 10 1 1_ 
LT napu тё 0.196 0.182. 0 195 0 196 0.188 0.194 0 189 0.199  (.I96 — Q.t03 
2.4 ойе» 224 — 0.162. D 165 0164 0.161  0O.t64 0 165 0.168 0.167 0.64 
3. M . nusihyerus 206 185 = 0.078 0077 0.073 0.073 0072 uma 01074 0.074 
3. M anhuiensis = 1 222 188 Ro = о 0001 0.069 0.072 0.068  U,D76 QUT? 0.074 
5 М anhurensa -2 223 187 88 L = 0.068 0.071 0.067 0.075 076 — 0.071 
6 M chen ugaster 4 1&4 83 79 TR = 0.071 0024 0.04  Q.U45 0.039 
7 M fiche 2a] ly? 33 32 Bl 24 = 0.015 0.044 бду 4.034 
8. M . lewogaster 1 188 82 17 To 27 7 == 0,041 пыл — 9.035 
9, M. b. byrngensi — E 227 192 83 87 86 50 5р 47 = aima ПОП 
10. M. b. byangensis — 2 224 190 84 88 87 51 51 48 5 == 0.011 
11. M. b. сарала 220 IB? 84 & a3 45 45 40 It 2 Еа 


PD Д аот оаа ЕО БЕСИН соз чш сш с ———= 


‘The nunibers helow the diagonal represent sequence substitutions, those above the diagonal represent ratio uf seguente divergem r 


11111111111122222222222333333333333444444444 
122344456669912224566789901223445788012223455999233446 778 
5218625870396916791025495846257398625421373518367958142172 


M.moschiferus CCCTAACTCCCCTCACCTACTCCTTTCCCTTCTACGTTCCCCACCCTCCACCTCCTC 
M.anhuiensis-1 T.T..GT...... T...CGTCTTCCCT.T.CTC.T..C.TTTGTT..TT ..CT.AT 
M.anhuiensis-2 T.T:aOT. 4e T...CGTCTTCCCT.T.CTC....C.TTTGTT. TT... CT. AT 
M.chrysogaster STTOGGICTIT.CT.;AC.T.:T-CC.. T STC GAs silo жеш TTCTA.. 
M.fuscus SOTT.OGG.ITI..CT.DPACOTOIOT.CC.STCSTG.s dem н ён GTTCTA T 
M.leucogaster STTC.IG...T..CT.TAC.TC.T.CC..TC. TC... Los mm C...TTCTA T. 
M.b.bijiangensis-i .TT.GG.. TT. T. TGTAC. T. TT, C. . T. C. TCG. . C. T.. TG. . T.. .TTCTT T. 
M.b.bijiangensis-2 .TT.GG..TT.T..GTAC.T..T.C..TTC.TCG..C. T.. TG. . T.. .TTCTT T 
M.b.caobangis .TT.GG..TT...T.TAC.T..T.C..TTC.TCG..C. T. . TG. .... .TTCTT.TA 
445555555555566666666677177701111117171171178888888888988808888259 
890244556788801134678900122234445677899011123455778888850C 
6218392879289378520544092129514565149251069270584902355]107 
M.moschiferus GTCCTTCATCGTCTATCCTCTATTAGCGTCCCTCTATATTCGCTCTTTCCCAGTTCTT 
M.anhulensis-1 A.TT.CT.C.ACTCGC..C.CG. ТА...... C....CTATCTC.CTTT.ACCACC 
M.anhulensis-2 A.TT.CT.C.ACTCGC..C.CG..,.TA...... C. ..CTATCTC.CTTT.ACCACC 
M.chrysogaster A.TTCCT.CT.C.C.C..C.CG.CGA. ACTTTCTCG.GC..A.C.CCC. T. C T.. 
M.fuscus A..T.CT.CT.C.C.C.TC.CG..GA.ACTTTCTC..GC..ATC, , CC. , S.C 
M.leucogsster A...CCT.CT.C.C.C. TC. CG. . GA. ACTTTCTC...C..ATC..CC. T .C. .. 
M.b.blijiangensis-1 AC..CC..CT.C.C.CT.C.CGA.G, TAC. TICTC.C.A..A.C..CC T ...TC. 
M.b.bijiangenais-2 AC..CC..CT.C.C.CT.CTCGA.G. TAC. TTCTC.C, A. . A. C. .CC.. T ... TC. 
M.b.caobangis AC..CC.GCT.C.C.CT.C.CGA.G, TAC. TTCTC. C. A. . A. C. .CC.. T... TC. 
111111111111111111111 
9999999999999999999000000000000011111111 
0112223566678899999002456667778800011222 
9254670703950134679589402580470614746259 
M.moschiferus CACCCACCTCCCTTTTATGTATCCACTTCTTICGACGCCTT 
M.anhulensis-1 TGTTT.TTCT.TC.CC..A...... ТС e G.ATT.C 
M.anhuiensis-2 TGTTT.TTCT.TC.CC. А...... TOi piae G.ATT.C 
M.chrysogaster T..T.G.TCT.T..CC.CACGCTT.TC.... TA...TT.. 
M. fuscus ...T.GTTC..T.CCC.CAC.CTTGTC. .CCTA.T.TT.. 
M.leucogaster T..T.GTTC..T.CCCGCAC.CTT.TC. .C. TA. .ATT., 
M.b.bijiangensis-l1 T..T..TT..TT..... САС.СТТ.ТССТ....... ттс. 
M.b.bijiangenais-2 Te :T. Ts ET ales, ыч САС.СТТ.ТССТ....... TTCC 
M.b.caobangis T..T..TT..TT..C..CAC.CTT.TCC. ,.. A. . ATTC. 


Fig. The variant sites (156 bp} of the eytocbrome h gene sequences ( 1 140 bp) of six species in musk deer 
The highlighted sites are those unique vunatlons in Anbu musk deer. The numbers abuse the sequences are the positions of those vartant 


sites 
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M.moschiferus 
M.chrysogaster 
M.fuscus 


M leucogaster 


M.b.caobangis 
M.anhuiensis-1 
M.anhuiensis-2 


A.alces 





Т.пари 


Fig.2 
The most parsimonious tree was constructed by using PAUP 3.0. The bootstrap values are shown above the branches (tree length = 502, СІ = 
0.787). The number vf one-step less parsimonious tree is zero. 


deer and the others is up to 7.3% on average. This 
value is equivalent to the largest sequence divergence 
(7.41%, Su et al., 1999) in musk deer and also 
falls in the spectrum of species sequence divergence of 
deer (4% — 12%) (Cronin, 1991). Hence. this te- 
sult, if confirmed at nuclear loci and supported bv eco- 
behavioral data, implies that the Anhui musk deer 
should be elevated to full species rank rather than a 
subspecies of either M . moschiferus ог M . berezovskii as 
suggested before (Wang et al., 1982; Wang et al., 
1993; Groves et al., 1995). According to the tree, 
Anhui musk deer branched out very early in Genus 
Moschus, which happened around 700 000 years ago if 
a silent divergence rate of 10% per million years in 
mammals were used {Irwin et al., 1991). 

[t is of interest to compare the molecular data with 
morphological data. Apparently, the molecular data 
revealed higher divergence of Anhui musk deer than the 
morphologic studies. at least on the maternal side. 
However, if we retrieve previous morphological stud- 


ies, Ánhui musk deer did show some morphological 
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M.b.bijiangensis-l 
M b.bijiangensis-2 
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The updated phylogenetic tree for Genus Moschus with Anhui musk deer added 


characteristic as described by Groves et af. (1995). 
"upper side with diffuse yellow spots in adult; rump 
without blackish bottom zone". Furthermore, in the 
multivariate analysis for measurements of skull and den- 
tal. one skull of Anhui musk deer was unexpectedly 
to fall outside the of both 


M . berezovskii and М. moschiferus in the discriminant 


found just range 
function plots {Groves et af., 1995). a morphological 
implication of its uniqueness . 

During the publication of this study. Li et al. 
(1999) reported their study on Anhui musk deer by se- 
quencing a 367 bp fragment of cytochrome b gene. 
They found 6.24% - 7.90% sequence difference be- 
tween Anhui musk deer and other musk deer species, 
which is consistent with our complete cytochrome b 


gene sequence data. 


Acknowledgements We are grateful to GOU 
Shi-kang, YU Yi-Chuan, ZHU Chun-Ling and QI Jin 


for their technical assistance . 


0000 http://Wwww.cqvip.com 


3 5 йй AS, RRMA DNA AHA b AED IR EFE ARES (ЖУ) 173 











References 


Cronin А О M, 1991. Mitochondrial DNA phylogeny of deer { Cervidae ) 
[J] J. Mamm. .72(37) 1553 - 566. 

Ellerman J R . Morrison-Scott T € 5. 1950. Checklist of Palaearclic and In- 
dian Mammals 1758 - 1946 [M / . Tonbridge. Kent, UK: Tonbridge 
Printers LTD 353 - 354. 

Felsenstein f, [985. Confidence шои» on phylogenies; an approach nang 
the boolstrap[ 1] . Evolution , 39:783 - 791. 

Felsenstein J, 1993. PHYLIP ; phylogeny inference package [ СР]. Version 
3, Se, Vinversily of Washington Seattle, USA. 

Gao Y T. 1963. Taxonomy of Chinese musk deer[],. Chinese J. Zool., 

15.479 - 487. 

Cau Y T, 1985 Classification and distribution of ihe musk deer 
í Moschus ) in Chuna, 4_ . In: Kawamichi T.Conlemporary Mammalo- 
gr on China and Japan "M. Tokyo: The Mammalogical Somely of 
Japan. 1143 ~ 1146. 

Groves С P, Wang Y X. Grubb P, 1995. Taxonomy of mush deer, genus 
moschus { Moschidae. Mammalia }! J]. Acta Theriologica. Sinica , 15 : 
181 ~ 197 

Groves © P, Feng Z J, 1986. The status of musk-deer from Anhui 
province ,China_J_ Acta Theriologica Sinia „61101 — 106. 

Grubb Г, 1982, The systemaues of Sino-Hymalavan musk deer 
(Muschus) , with particular reference to the species desenhed by 
Hodgson BH[J_ . Saeugetierkundliche Mitteilungen ,30:127 — 135 

lean D M, Kocher T D, Wilson A C, 1991. Evolution of the cytuchrome b 
gene of mammals[]].J. Mol. Evol ..32: 128 - 144. 

Komar S, Tamura K, Ne М, 1993. MEGA; molecular evolutionary genetics 

analysis! CP]. Version 1.0. The Pennsylvania State Umversity, Uni- 


versity Park PA 16802. 

Li M Li Y G, Sheng Н L.1999. Taxonomy status of Апаи musk deer re- 
visited[]]. Chinese Science Bulletin .44(2); 188 ~ 192 

la Z X,1981 А new species of Chinese musk deer J] Zidogucal Res , 
2:157 ~ 161 

«шаб N, Gao Y T. 1990. А review of the distribution af all species uf 
deer { Tragulidae, Moschidae and Cervidae ! зд China: J|. Mammal 

Review , 20:125 ~ 144. 

Sheng Н L.1989 Moschidae (Chapter 5)[ A]. In: Sheng H L. The Deer 
m China[ M] Shanghai; Eastern China Normal University. Press. 45 
~ 88. 

Su В , Wang Y X, Lan Н et af, 1999. Phylogenetic study of complete 1y- 
tochrome H genes of mush deer {Genus Moschus using museum 
samples( ]] . Molecular Phylogenetics and Evolution ‚12241 ~ 249. 

Swofford D L, 1989. PAUP; phylogenetic analys using parsimony | CP ] А 

Version 3.0.lllinosis Natural History Survey. 

Walsh P S. Metzger D А, Higuchi R,1991. Chelex 100 a. a medium tor 
simple extraction of DNA for PCR-hased lyping frum forensic maten- 
al J". Brot. Techniques ,10:506 - 513. 

Wang Q S. Hu X 1. Yan Y Н. 1982. А new subspecies of Moschie 
moschiferus эп China; М. m. anbaiensis [I]. Acta ТАетийо аса Sint- 
va .2(2):133 - 138. 

Wang Y X. Ma S L, Lı C Y.1993. The taxonümy . distnbulion and status of 
furest musk deer in China[ A ). 1n; Ohtaishi N,Sheng H L Deer of 
China;Biology and Management [M |. Tokvo: Elsevier. Tokyo. 22 — 
30 


SC ОМА ARE XC b BAS KASD HET 


a Edu 
(ТЕЕ ре BA 20 ЗЕ АТ EH 650223) 


ЕЕ #7 


жш? 
OTRE „е ХЕК BE A Ж fe РЕ Г.) 


СКЕЧ TREE RHE 230039) 


HE. BH 8O4E(RR ЛЕ [Moschus ( moschifer- 
us / berezovskii ) anhuiensis | Ж УНЫ Ж, Bob EE 
іу — Efi. EAE. RU SCENES 
БЕК ait IT SAK DNA MES b ASK 
PIs, ВРА НН, RRA RR PH he 
Bst fe 2 ILC. ES THO BE BU DT GRIEF PMR 


Hj. RRR PHAR. TER] LC f Rp Н) 
DNA ЕЗД ЗЕ SE CARROL, ВК, E 
ЖЖ DNA FE SUR SESS Sc SENE SE MEN ALS TR — 
APP (Moschus anhuiensis), M EB AAA 
Бану oO ЧЕ Ж 


Xu. B BR: LRM, #НИ@ ЖЪН, TRAE 


rH. 0959.842, 0349-.5 ЖЕЕП. А 


ХШ =: 0254 - 5853( 2001 103 - 0169 - 05 


